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It is well known that quasi-local black hole horizons depend on the choice of a time coordinate 
in a spacetime. This has implications for notions such as the surface of the black hole and also on 
quasi-local physical quantities such as horizon measures of mass and angular momentum. In this 
paper, we compare different horizons on non-spherically symmetric slicings of Vaidya spacetimes. 
The spacetimes we investigate include both accreting and evaporating black holes. For some simple 
choices of the Vaidya mass function function corresponding to collapse of a hollow shell, we compare 
the area for the numerically found axisymmetric trapping horizons with the area of the spherically 
symmetric trapping horizon and event horizon. We find that as expected, both the location and 
area are dependent on the choice of foliation. However, the area variation is not large, of order 
0.035% for a slowly evolving horizon with m = 0.02. We also calculate analytically the difference in 
area between the spherically symmetric quasi-local horizon and event horizon for a slowly accreting 
black hole. We find that the difference can be many orders of magnitude larger than the Planck 
area for sufficiently large black holes. 
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cr. 

^ . 



> 

O 
0^ 
0^ 
(N 

l> 

O 
O 



X 



I. INTRODUCTION 

An important feature of black hole spacetimes is the 
existence of trapped and marginally trapped surfaces. 
These surfaces are a key element in the Penrose singular- 
ity theorem and are frequently used in numerical relativ- 
ity to indicate the existence of a black hole. Recent work 
has introduced related concepts such as trapping horizons 
and dynamical horizons both of which, ignoring certain 
technical conditions, can be viewed 3-surfaces obtained 
by time-evolutions of marginally trapped surfaces. It has 
been shown that these surfaces have many properties sim- 
ilar to event horizons, in particular their thermodynamic 
properties [1-3]. Trapping and dynamical horizons can 
also be used to extract data about the black hole, such as 
its mass and angular momentum, directly from the strong 
gravity region without needing to rely on the asymptotic 
behaviour or exact isometrics [1]. 

However, it has been known for a long time that the 
location and existence of trapped surfaces depends on 
the choice of spacetime foliation, i.e. on the choice of 
the time coordinate, and this is the main reason why 
event horizons have often been preferred over trapping 
and dynamical horizons as a definition of the surface of a 
black hole. In particular this means that any given, suit- 
ably general, spacetime that contains a trapping horizon 
is likely to contain, in principle, an infinite number of 
trapping and dynamical horizons that intersect but are 
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distinct [5]. The key results for quasi- local horizons, for 
example the laws of black hole mechanics and the energy 
and angular momentum balance laws [6] , are sufficiently 
general and are applicable to all of these different hori- 
zons. Nevertheless, if we wish to assign say mass and spin 
to the black hole quasi- locally, the time evolution of these 
quantities will depend on which horizon we choose to use. 
Thus, to use quasi-local horizons to extract meaningful 
data from black hole spacetimes it is important to in- 
vestigate the location of these distinct horizons and how 
this affects their related properties such as mass, angular 
momentum, linear momentum etc. 

This issue is potentially relevant for gravitational wave 
astronomy when one seeks to infer the physical parame- 
ters of the source from the observed gravitational wave 
signal. For example, in modeling binary black hole coa- 
lescence it becomes important to combine results from 
post-Newtonian theory and numerical relativity. Nu- 
merical relativity can in principle evolve a binary black 
hole system starting from the inspiral phase all the way 
through the coalescence (see e.g. [7] for a recent review). 
However, this becomes computationally impossible if we 
start the black holes very far apart, and if we wish to 
cover a large parameter space. Post-Newtonian theory 
on the other hand treats the black holes as point particles 
endowed with a mass and spin, and solves the field equa- 
tions in powers of v/c with v being the orbital velocity of 
the black holes. This approximation works best when the 
black holes are far apart and breaks down near the coales- 
cence. Thus, in order to fully solve the binary black hole 
problem it becomes essential to combine the two frame- 
works. The literature on this topic is large and a full 
discussion is beyond the scope of this paper; see e.g. [8- 
14]. For our purposes, we note that the two frameworks 
treat black holes very differently; post-Newtonian theory 



treats the black holes as point particles while numerical 
relativity deals with the horizons non-perturbatively. It 
turns out that using different flavors of Post-Newtonian 
approximants leads to biases in the values of physical 
parameters [15] (as infered from the gravitational wave- 
forms). Numerical relativity on the other hand deals with 
black hole horizons, and if we decide to use quasi-local 
horizon measures for calculating physical quantities, the 
values of the black hole parameters will be affected by the 
choice of foliation. This effect is absent when the black 
holes are isolated; thus we expect it to be negligible in 
the inspiral phase and to be larger in the dynamical co- 
alescence phase. In comparing the results from the two 
frameworks, it is useful then to quantify this source of un- 
certainty and to show that errors in the physical param- 
eters are smaller than the biases between, say, different 
post-Newtonian approximants. 

Another area where the location of horizons is impor- 
tant is in black hole thermodynamics. If the area of 
a black hole horizon is to play a role analogous to en- 
tropy via the Bekenstein-Hawking relation it is impor- 
tant to know which area this relation should be applied 
to in which situations. For stationary black holes, cross- 
sections of the event horizon are marginally trapped sur- 
faces and thus the different notions of black hole hori- 
zons (i.e. event horizons and trapping horizons) coincide. 
However, in a dynamical spacetime the event horizon will 
not coincide with the location of any of the trapping hori- 
zons. In general they will have different areas on a given 
spacelike hypersurface [16, 17]. Furthermore, although 
the location of the event horizon does not depend on the 
choice of spacetime foliation, when the event horizon is 
growing or shrinking, its area will depend on the choice 
of foliation. 

Non-stationary spacetimes with dynamical black holes 
are complicated, and only a few exact dynamical black 
hole solutions are known. One class of solutions is the 
Vaidya solutions, which describe the evolution of a spher- 
ically symmetric, radially moving, pressureless null fluid 
with a freely speciflable mass function m(v). This is a 
useful and non-trivial toy-model for a dynamical black 
hole. These spacetimes can be used to qualitatively 
model the physical accretion of matter by a black hole. 
Some authors [18, 19] have even suggested using the 
Vaidya spacetimes to simulate the decrease in area of 
a black hole due to Hawking radiation, by allowing the 
infalling matter to be negative. 

The uniqueness of dynamical horizons in general space- 
times was investigated in [5] where it was shown that 
the foliation of a given dynamical horizon by marginally 
trapped surfaces is unique, and that it is not possible to 
foliate a region with dynamical horizons since they inter- 
sect. In the Vaidya solution it has already been shown ex- 
plicitly that marginally outer trapped surfaces (MOTS) 
can be found at various different locations depending on 
the slicing [20], although this reference did not investi- 
gate how the parameters of the black hole, such as area, 
mass and angular momentum vary with slicing and did 



not consider the case of timelike trapping horizons. Ben- 
Dov [21] showed how to find marginally outer trapped 
surfaces arbitrarily close to the event horizon. It was 
also shown that there arc flat regions inside the event 
horizon (the region contained within the shell) where no 
trapped surface passes. Bengtsson and Senovilla [22] an- 
alytically constructed trapped surfaces that pass through 
the flat region. This is not a contradiction; while there 
are some flat regions where no trapped surface can be 
located, there are trapped surfaces in other parts of the 
flat region. In spherically symmetric spacetimes spheri- 
cally symmetric marginally trapped surfaces are easy to 
find and in the case of Vaidya are given by the condition 
r = 2m{v), where r is the areal coordinate (r = y^A/Wn 
with A being the area of the MOTS) of the surfaces of 
spherical isometry. It is also easy to show that such sur- 
faces foliate a trapping horizon which we will refer to 
as the spherically symmetric trapping horizon. Due to 
results in [5], any other dynamical horizon will lie par- 
tially outside the spherically symmetric dynamical hori- 
zon somewhere. For example, in [23] analytic solutions 
were presented where closed trapped surfaces extend into 
the region between the spherically symmetric trapping 
horizon and the event horizon. The surfaces we find here 
do the same. 

In this article we will look at the variation of black 
hole parameters for horizons located on various different 
slicings of the same spacetime. We will examine several 
different spherically symmetric mass functions; a linear 
mass function designed to see the behaviour for slowing 
evolving black holes and a tanh log mass function, de- 
signed to see the behaviour in a short collapse of a shell to 
a black hole. We compare marginally outer trapped sur- 
faces found on non-spherically symmetric hypersurfaces 
with those found on the spherically symmetric slicings. 
We are able to locate the spherically symmetric horizons 
analytically but the non-spherically symmetric horizons 
are located numerically using the AHFinderDirect thorn 
[24, 25] of the Cactus framework [26, 27]. A number 
of useful analytic relations for the Vaidya spacetime are 
given in the appendix. Background detail on the proper- 
ties of trapping horizons and their thermodynamics can 
be found in [1] and [3]. 



II. THE VAIDYA METRIC 

To examine the slicing dependence of quasi-local hori- 
zons we need to consider dynamical spacetimes. The 
Vaidya solutions have a number of nice properties that 
make them popular for investigations of this type [20-22] . 
The Vaidya spacetimes [4] are a class of spherically sym- 
metric, non-vacuum spacetimes with line element in ad- 
vanced null Eddington-Finkelstein coordinates (u, r, 9, (p): 
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The function m{v) is a a freely specifiable mass func- 
tion; it coincides with the Misner-Sharp mass [29] in 
this case. In advanced Eddington-Finkelstein coordi- 
nates the metric is well defined across future horizons. 
The Vaidya metric solves the Einstein equations with an 
energy-momentum tensor of the form 



Tab = Enam, 



(2) 



where n" is an ingoing radial null direction and i<^ is a 
function that depends on the normalization of n". For 
example, for the case where Ua = —daV then we have 
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The Vaidya so- 



lution can be interpreted as describing the radial collapse 
of pressureless null dust (i.e. infalling radiation). From 
the appendix we can see that, in a spherically symmet- 
ric foliation, trapping horizons will occur at r = 2m{v) 
provided we keep m{v) > 0. 

Here we will study various mass functions for m{v). 
The simplest is a linear mass function of the form 



m{v) ~ rriQ + mv. 



(4) 



The linear mass function is suited to situations where the 
black hole is accreting matter at a constant rate (for a 
finite time duration to ensure that the mass stays finite). 
While the assumption of spherical symmetry is some- 
what artificial, real astrophysical black holes do have very 
small mass accretion rates whether from the accretion of 
surrounding stars and gas or purely from the accretion 
of photons from the cosmic microwave background. In 
the Vaidya solutions the Misner-Sharp mass on each in- 
going constant v surface is a constant. This refiects the 
fact that the mass is flowing inwards at the speed of light 
and the mass contained within a shell of radius r is con- 
stant as the radius decreases. Although the spacetime is 
dynamical we can still define observers who remain at a 
fixed areal radius r and fixed 9 and (f). In terms of the 
proper time r of such observers we have 
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■\/l — 2m{v)/r 



(5) 



At large distances from the black hole r ^ m{v) we can 
use this to relate directly the mass flux in terms of the null 
coordinate v to the mass flux as seen by constant r, 9, (p 
observers who would be static observers in an exactly 
static spacetime. 



dm{v) dm^v) 



dv 



dr 



(6) 



We can give a very rough indication of the order of mag- 
nitude for m that might be expected for certain astro- 
physical cases. The Eddington rate is used to estimate 
the maximal rate at which infalling matter can be sup- 
ported by its own radiation pressure and large, luminous 



black holes are typically found with luminosities between 
10 and 100% of the Eddington limit [30]. For a black hole 
accreting at a tenth of the Eddington limit the dimen- 
sionless accretion rate is approximately [16] 
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This is ^ 10"^* ergs per second for a solar mass black 
hole. The matter falling into the black hole is ten times 
the energy being emitted as light. Since this accretion is 
usually associated with a disk though it will not be spher- 
ically symmetric. For a black hole accreting purely from 
the Cosmic Microwave Background, which is assumed to 
be isotropic, we have from the Stefan-Boltzmann law, 
approximately 
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(8) 



For a black hole whose dynamics are dominated by evap- 
oration through Hawking radiation we have 
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(9) 



For numerical purposes we will investigate evolutions 
with mass rates much higher than astrophysical rates, 
typically jmj ^ 0.01. To the extent that 2\/to <t^ 1, 
the spherically symmetric trapping horizon will still be 
a slowly evolving horizon in the sense of [32]. We will 
also look at mass functions of the tanh log form, with 
m{v) — for w < and for w > 



m{v) 



mo 
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1 + tanh I log 
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T^ 



(10) 



The first derivative of this tanh log mass function van- 
ishes at w = 0. Therefore the metric and its first deriva- 
tive are contiunous at w = 0. For T and Too greater than 
zero, this mass function models the collapse of a hollow 
spherical shell of matter m{v = 0) = that asymptoti- 
cally settles down to an isolated black hole of mass to-q for 
w — > oo. These types of mass functions model situations 
where the black hole grows initially quite rapidly and 
then asymptotically settles down to its final static state. 
The mass function reaches half its asymptotic value when 
v = T. 

The maximum of m for these tanh log functions occurs 
at w = T/v3 and takes a value ^ mo/T. Therefore these 
solutions will not be slowly evolving in the sense of [32] 
for TTiQ and T of similar sizes. 



A. Location of spherically symmetric horizons 

The trapping horizons are three-dimensional surfaces 
H, foliated by closed spacelike two surfaces for which the 
future directed null normals £"- and n" satisfy 
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Cn9(^i) <0. 



(11) 



Here 0(£) and 0(„) are the expansions of £" and n" respec- 
tively, and C\ is the Lie derivative along n". Dynami- 
cal horizons are also three-dimensional surfaces foliated, 
as above, by spheres with 0(£') = 0, 6'(„') < 0. How- 
ever, CnOit) < is replaced by the requirement that H 
be spacelike. For the spherically symmetric horizons in 
Vaidya, these notions coincide [1]. 

For spherically symmetric slicings, the null normals as- 
sociated with the trapping horizon will be radial null vec- 
tors and the location of the trapping horizon is just given 

by 



r ~ 2m{v) . 



(12) 



This is a spacelike surface for m > 0, a null surface for 
TO = and a timelike surface for m < 0. For the linear 
mass function, in terms of the timelike coordinate t — 
V ~ r, the horizons will be located at 



2to,o + 2mt 
I -2m 



(13) 



For the tanh log mass function the horizons are located 
at the solution of the cubic 

r^ + (2i - 2too) r^ + (t^ + T^ - 4moi) r - 2TOot^ = 0. 

(14) 
The equation for the horizon r — 2m{v) has a single 
unique solution on each surface of constant v. But it 
can have multiple solutions, corresponding to multiple 
horizons, on surfaces of constant t. For to-o > this 
cubic function always has at least one positive real root 
and guarantees that there will always be at least one 
horizon. In fact, it can be shown that there will only be 
a single horizon for t > 2to-o — ■\/4to,q — T^. For small 
values of t there are multiple horizons and this will also 
be the case for other mass functions. 

We will use numerical methods to solve the trapping 
horizon equations for non-spherically symmetric slicings 
and postpone discussion of them until the next section. 
Due to results in [5] we expect that the non-spherically 
symmetric trapping horizons will intersect the spherically 
symmetric ones. 

The event horizons are defined as the past causal 
boundary of future null infinity and are generated by null 
geodesies that fail to reach infinity. The event horizon is 
always a null surface since it is a causal boundary. In the 
Vaidya spacetimes it is generated by radial outgoing null 
vectors that satisfy 



dr 1 / 2m(v) 
dv ^ 2\ r 



(15) 



This first order ordinary differential equation generates 
the path of all outgoing radial null geodesies. In order to 
give the location of the event horizon it requires a bound- 
ary condition that corresponds to the known location of 
the event horizon at some particular point. In practice 
this is usually given by the position of the event horizon 
at some future point, either when the black hole evapo- 
rates entirely or settles down to a stationary state. If the 



black hole at some point settles down to a Schwarzschild 
black hole with no further matter accreting, then the 
event horizon can be located by tracing back the null 
rays from the future Schwarzschild radius. However, in 
the situation where the black hole is accreting matter at 
a steady rate and is a suitably long way from changing 
to a different state one can find the approximate location 
of the event horizon by imposing the condition 
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(16) 



This just refiects the fact that the event horizon is grow- 
ing at a steady rate [16]. In this case, equation (15) has 
the general solution 



m{v) 

f = ; — 

4to 

For m <C 1 this gives 



1 — Vl — 16to 



r = 2m{v) (l + 4to + ©(to^)) 



(17) 



(18) 



and thus we expect that the event horizon will be outside 
the spherically symmetric trapping horizon for m > but 
inside for m < 0. For a solar-mass black hole accreting 
at a tenth of the Eddington rate the difference in areas 
between the event horizon and the spherically symmetric 
trapping horizon will be around 10^^ in units of Planck 
area, while for a supermassive black hole of mass 10^ solar 
masses, accreting purely form the CMB, the difference in 
areas will be around 10^" in Planck units [16]. In terms 
of the coordinate t = v — r for the linear mass function 
the event horizon has radial coordinate 
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1 -2to 



(19) 



This is just the location of the spherically symmetric 
trapping horizon with a constant offset of 8toto provided 
TO ^ 1 . In this approximation the generators of both the 
trapping horizon and the event horizon have the same 
components but the norm of the generators is 4to for the 
trapping horizon and zero for the event horizon. 



B. An axisymmetric spacetime slicing 



We consider a simple axisymmetric slicing of the form 



t 



(20) 



where z — r cos 9 and a is a parameter that determines 
how far away from spherically symmetry the constant i 
surface is. We reserve the symbol t for spherically sym- 
metric surfaces where a = 0, which just gives the usual 
Eddington-Finkelstein time coordinate (although not the 
Schwarzschild time coordinate). Hypersurfaces of con- 
stant t are always timelike for |a| < 1 since the normal 
to a given hypersurface, i", has norm 
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2TO(f ) 



(l + a cos 9 + 0^ cos^ 



(21) 



On the slice i = we have w > everywhere for |a| < 
1. In addition, since (l + acos0 + a^ cos^ 6*) is always 
greater than 3/4, for real a this hypersurface will become 
timelikc near the horizon for a > -\/7/2 ^ 1.323. 

On each slice with a constant value of i, a two dimen- 
sional marginally outer trapped surface (MOTS) can be 
searched for, satisfying just 9ti\ = 0. In some spacetimes 
there may be multiple MOTS on a given hypersurface. 
This will not occur in the Vaidya spacetimes with linear 
mass function since the horizon condition, r — 2m{v), 
is linear in r. The marginally outer trapped surfaces 
can typically be stacked to form three dimensional trap- 
ping horizons, provided the other conditions 0/^) < and 
"°Va^(„) — are satisfied too (or a dynamical horizon if 
the resulting surface is spacelike). 

For the case a = the slicings will be spherically sym- 
metric. The orbit of spherical rotations of each point will 
lie entirely in the hypersurface. This will lead to MOTS 
that are spherically symmetric, where each surface is just 
the orbit of spherical rotations. Commonly used slicings 
in numerical relativity such as constant mean curvature 
surfaces will typically be spherically symmetric in the 
Vaidya spacetime. In this case the horizons will be lo- 
cated uniquely at r = 2m(v) and every point on the 
surface will have the same value of the advanced time v. 
This will not be the case for surfaces with a 7^ 0. In these 
cases the mass parameter m{v) will not necessarily take 
the same value on different points of the surface. The 
MOTS will then extend into regions where the matter is 
more compact and regions where it is less compact. Sur- 
faces with a > will be the same as surfaces with a < 0, 
with the north pole {6 = 0) interchanged with the south 
pole {6 — tt) or, equivalently cos9 — > cos(7r — 9). 

Each different choice of i will lead to different MOTS 
than can be stacked to form different trapping horizons. 
Each of these trapping horizons will have an associated 
natural foliation, the one that makes each two dimen- 
sional surface a MOTS. However we can also consider the 
intersection of a given trapping horizon with a spacelikc 
hypersurface of a given a value different from the natural 
one. We can compute the area of the these closed two 
spheres, although for a 7^ they are not MOTS. For 
the spherically symmetric trapping horizon (SSTH) in 
Vaidya spacetime with linear mass function and m <C 0, 
the area will be 

^SSTH = 47r (2mf ( 1 + Arfi + 2m— + -ma^ ) . (22) 
V m 3 / 



spherically symmetric surface with a = but the same 
constant value of t we see that the area will be greater if 
TO > and smaller if m < 0. 



C. Horizon deformations 

The choice of different foliations leads to different 
MOTS that can be deformed into each other. Using 
these deformations it is possible to examine how certain 
geometric properties of surfaces change with the defor- 
mation [31]. For example, it is possible to see how the 
area element, e, changes as MOTS are deformed into one 
another. 

We can search for MOTS on slices with different val- 
ues of t but because these marginally trapped surfaces do 
not lie on the same hypersurface it is difficult to compare 
them. If they do not intersect (and typically they don't) 
it may or may not be possible to find a trapping horizon 
that connects both surfaces. We want to distinguish be- 
tween marginally trapped surfaces that are evolutions of 
one another and marginally trapped surfaces that cannot 
be evolved into one another along a trapping horizon. 

For a given choice of t but two different values of a the 
two slices thus defined will intersect when z = 0. The 
slice i = with a = does not lie entirely to the past 
or future of the slice i = with a = 1. However, the 
horizons on these two slices will typically not intersect. 
On the two dimensional surface where these two slices 
do intersect, z = 0, the horizon for a = is found to lie 
outside the horizon with a 7^ 0. 

We can define a vector field, X, on surfaces that can 
be used to define variations of this surface, 8x ■ Following 
[32] we can write this vector field on a marginally trapped 
surface in terms of the two null normals to the surface as 



X" = Sr - Cv!^^ 



(24) 



where B and C are functions on the surface. For deforma- 
tions of the spherically symmetric surface where A = 
we have 



--4-4 



(25) 



for the Eddington-Finkelstein coordinates t and r. The 
variation of intrinsic geometrical scalar properties can be 
calculated using 



A similar calculation can be performed for the intersec- 
tions of the event horizon with the constant t surfaces. 
Taking the position of the event horizon to be (18) and 
again the approximation to <C wc find 

Aeh = 47r (2to)^ ( 1 -h 12to -H 2to V -mo? ) . (23) 

\ TO, 3 / 

If we compare the area of the event horizon on a constant 
t hypersurface with the area of the event horizon on a 



^x^ = C.x4>, 



(26) 



and for general tensors, by projecting the Lie derivative 
onto the surface 

^XWab == qlq^CxWab- (27) 

The variation of the area element e for example satisfies 



<5xe = -C9(n-^e 



(28) 



Variations of extrinsic properties are a little harder to 
calculate. We are particularly interested in deformations 
of the form 



^X»(l) 



0, 



(29) 



since this will generate a class of marginally outer 
trapped surfaces. The calculations in [32] give 
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B<5(,)0(,), (30) 



where d is the covariant derivative compatible with the 
intrinsic metric of the two-surface. Notice that for radial 
f " we have (3) and thus 



5(1)^ it) = 
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since the shear is zero. Furthermore, axisymmetric vari- 
ations from the spherically symmetric surface will satisfy 
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C - 2mB = 0. 



(32) 



In the case where both B and C are constant we get 
the evolution along the spherically symmetric trapping 
horizon. In [32] it was assumed that the deformations 
were I -oriented such that B > Q everywhere. But another 
solution of this equation is 



B 



2m 



C = —kcosB, 



(33) 



(34) 



where /c is a constant of integration. Since the variation in 
this case is dependent on 9 it will not generate spherically 
symmetric marginal surfaces, and since B = C = for 
6 = tt/2 this will not generate evolution along a trapping 
horizon. The area change in this case is easily computed 
as 
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(n)« 



0. (35) 



In this case, at least, the area of the spherically sym- 
metric surface is extremal. For a general axisymmetric 
variation expanded as a Fourier series, the area of the 
spherically symmetric surface is extremal provided the 
coefficients a/^) of the cos(n^) terms in C are zero for 
even n. 



III. NUMERICAL RESULTS FOR 
AXISYMMETRIC SLICINGS 

We now compare the different trapping horizons one 
obtains for the axisymmetric slicings given by (20) where 
we consider the spherically symmetric mass functions (4) , 
(10). For the numerical calculation in the following we 
use, as in previous work [20], the Cactus computational 



toolkit [26, 27] which is widely used in numerical rela- 
tivity, for example to perform binary black hole, neutron 
star, or stellar core collapse simulations. For each leaf we 
write the 3+1-decomposed Vaidya metric onto the Carte- 
sian grid of Cactus. In each 'time step' we write the next 
leaf of the slicing onto the Cartesian grid. This mimics 
a numerical time evolution of an initial Cauchy slice of 
a dynamical BH spacetime. To locate the trapping hori- 
zon at each 'instant of time' we use the apparent horizon 
finder AHFinderDirect [24, 25]. This way we obtain the 
location and area of the trapping horizon as a function 
of time. In the following we compare these functions for 
the different trapping horizons we detect on the axisym- 
metric slicings from above. 



A. Location of tilted versus untilted horizons 

The MOTS lying on the hypersurfaces i = v — r — 
az = with different values of a do not intersect. An 
examination of the marginally outer trapped surface with 
a =- 0.83 shows that the points at the north pole = 0, 
equator 9 = -k 12 and south pole 9 = tt have the four- 
dimensional spacetime coordinates given in Table (I) 



TABLE I: Coordinate location of various points on surface 
with t = 0, a = 0.83 and m{v) = 1.0 + O.Olu 



surface point 


r 


V 


t 


zmax (north pole) 

xmax (equator) 

zmin (south pole) 


2.0454 
2.0386 
2.0247 


3.7498 
2.0386 
0.3374 


1.7044 
0.0000 
-1.6872 



Trapping horizons can be formed by stacking surfaces 
found on different hypersurfaces with different values of 
i but the same value of a. These trapping horizons, each 
with their own value of a will intersect one another. This 
can also been seen in Table (I). The points on the equator 
lie inside the horizon located on the spherically symmet- 
ric slicing t ~ v — r — Q, where every point on the horizon 
has r = 2.0408 and t = Q. The point at the north pole 
lies inside the spherically symmetric horizon on the slice 
t = 1.7044, which lies at r = 2.0756, and the point at the 
south pole lies outside the spherically symmetric horizon 
on the shce t = -1.6872, which lies at r = 2.0064. 

In terms of the deformations discussed above (25) we 
see that to deform the spherically symmetric surface at 
t = we need B > and C > at the north pole and 
B < and C < at the south pole, where in both cases 
the size of C is much smaller than B. 

Figure (1) shows the locations (r coordinate) of the 
north pole, equator and south poles for MOTS against a 
for a mass function of the form m{v) = 1.0 -I- O.Olw. For 
values of a less than zero the location of the MOTS is 
the same with the north pole 9 = and south pole 9 = tt 
interchanged. Near the spherically symmetric MOTS at 
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FIG. 1: The location (r-coordinate) of the north pole, equator 
and south pole versus diflferent values of a for a linear mass 
function m{v) = 1.0 + 0.01«. The positions for negative a are 
the same as for positive a with the north pole and south pole 
interchanged. Near a = there is no change in the position 
of the equator and the change in the position of the north 
pole and south pole are equal but opposite in sign. 



a — the change in the r coordinate of the equator 
is zero, whereas the change in the r coordinates of the 
north and south poles is equal but opposite in sign. This 
supports the conjecture that for these cases C is of the 
form (34). 



B. Area dependence on a 




FIG. 2: Slicing dependence of the area at i for the linear mass 
function (top) with m — 1.0 and m — 0.02 and for the tanh 
log mass function (bottom) with m = 1.0 and T = 1.0 



The results for the area dependence of a are shown 
in Fig. (2) for the linear mass case and Fig. (2) for the 
tanh log mass function. The a = surfaces have the 
largest area. As expected the area of the tilted horizons 
is unchanged under a change in sign of a. The difference 
between the area on a = shces and a 7^ shces depends 
on rh as shown in Fig. (3) for the linear mass function. It 
is zero for m — 0, where the spacetime just reduces to 
the Schwarzschild solution and the horizon is an isolated 
horizon. The difference between the area of tilted and 
untilted horizons increases with increasing magnitude of 
rh. The difference is slightly larger for positive values of 
TO than for negative values. This difference in area does 
not seem to follow a simple power law and it is difficult to 
extrapolate what the difference in areas will be for values 
of TO different form the values examined. However, the 
difference in areas for values of m much lower than the 
lowest positive value of to considered here (to — 0.02) 
are likely to be much less than the ^ 0.035% difference 
found here for that case. 

Finally, Fig. 5 compares the area of the event horizon 
(or more precisely, the intersection of the event horizon 
with the spatial slices corresponding to different values 
of a) with the area of the MOTS as functions of a. It is 
interesting that the variation in the event horizon area is 



much larger than for the apparent horizon. 



C. Dependence of rate of change of area on a 

For the linear mass function the a 7^ surfaces seem 
to have the same rate of change of area as the a = 
surfaces, at least to within numerical accuracy. For the 
tanh log mass function, as depicted in Fig. (4), the area of 
tilted horizons grows at a greater rate than the untilted 
areas and both converge to the asymptotic isolated hori- 
zon area, which will also be the asymptotic area of the 
event horizon. In this sense the spherically symmetric 
trapping horizon plays a special role as it is the one with 
the smallest rate of increase of area. This would seem 
to contradict the conjecture made in [33] that the pre- 
ferred trapping horizon should be the one where the area 
increases the most and approaches the event horizon the 
fastest. However, it is not known whether this behaviour 
is repeated in more general situations with other mass 
functions or other slicing conditions. 
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m 

FIG. 3: Area difference for two slices in different Vaidya space- 
times witli linear mass function and m = 1.0, The difference 
between the areas increases for large m. The untilted areas 
are larger than the tilted areas for m both positive and nega- 
tive. However, the difference is much larger for rh > 0. Only 
for rh = do the areas coincide. 
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FIG. 5: Area of the event horizon A^h (dashed line) and 
the MOTS Aah for different values of a for the linear mass 
function with rh = 0.02. The y-ajcis on the left refers to the 
area of the event horizon and the axis on the right to the 
MOTS; the scale of Aah has been expanded to better show 
its variation. The variation of Aeh is much larger than the 
variation of Aah and the variation of Aah is small relative to 
the difference between Aeh and Aah- 




Time 

FIG. 4: Area versus time t for a — and a = 0.83 for the 
tanhlog mass function with m = 1.0 and T — 1.0. Both areas 
converge rapidly to the asymptotic isolated horizon area IGvr. 
The tilted area grows faster since it starts from a lower point. 



IV. CONCLUSION AND DISCUSSION 

We have investigated the slicing dependence of the ap- 
parent horizons in the Vaidya spacetime. For a given slice 
we have looked at three types of horizons; the MOTS that 
lie on the slice, the intersection of the slice virith the event 
horizon and the intersection of the slice with the spher- 
ically symmetric trapping horizon. We have examined 
rather simple axisymmetric slicings in this simple space- 
time. These simple slicings do not exhaust all possible 
axisymmetric slicings and there are many other possible 
slicings that are not axisymmetric. 

ExpHcitly the location of the horizons is different for 
different choices of horizon and different choices of slicing. 
However, for slowly evolving horizons, the areas do not 



vary by much, although the larger the rate at which the 
black hole accretes matter the larger the difference in 
the areas. For the m — 0.02 Vaidya solution with linear 
mass function the area can vary by approximately 0.03% 
as a is varied from to 1. Beyond a = 1 the Cauchy 
surface becomes timelike for certain values of 9. This is 
the smallest value of m that we investigated but it is still 
much larger than the mass accretion rate expected for 
astrophysical black holes. 

From a purely practical point of view we've found ev- 
idence that the parameters of the black hole, derived 
from its trapping horizon, do not change significantly 
when looking at reasonably simple slicings. It is still un- 
known whether the properties would change drastically 
for certain unusual slicings, but these are unlikely to oc- 
cur in normal numerical simulations. It is also reassuring 
that these differences are much smaller than the param- 
eter biases of up to 10% found between different post- 
Newtonian models [15]. This variation of the parameters 
with foliation will also occur if one uses the event horizon. 

From a conceptual point of view, we have demon- 
strated explicitly the known result that the location of 
the black hole surface and some of its properties such 
as the area depend on the choice of the spacetime slic- 
ing. There are a variety of different possible responses 
to the issue of non-uniqueness of the trapping horizons 
in a given spacetime. The oldest approach is to focus 
purely on the event horizon as the unique indicator of 
the black hole and its properties. The location of the 
event horizon is independent of the choice of foliation. 
This approach however causes trouble in certain quan- 
tum inspired spacetimes [3, 18] where no event horizon 
exists. In numerical settings, the acausal definition of the 
event horizon means that it is only known after a simula- 
tion is finished, preventing it from being used to analyse 



the simulation's state as it proceeds. 

Another approach is to accept all horizons on an equal 
footing as purely a property of the geometry. In this pic- 
ture it is not clear how to associate unique properties to 
the black hole such as a horizon area or horizon angular 
momentum. Because the horizons intersect one cannot 
rely on specifying, for example, the outermost trapping 
horizon. However, it may be possible to formulate a gen- 
eralized second law for each possible foliation of space- 
time and in this context use the horizon defined by the 
chosen foliation. 

A third approach is to look for the boundary of the re- 
gion that admits trapped surfaces. This surface should be 
spherically symmetric in a spherically symmetric space- 
time. If this surface lies strictly outside the spherically 
symmetric trapping horizon then, by the results of [5] , it 
cannot itself admit the structure of a dynamical horizon. 
This surface will have a location that is independent of 
any given foliation, although it may not have a simple 
thermodynamic interpretation. 

A fourth approach is to look for properties that select 
out a certain preferred trapping horizon in the space- 
time [34], such as the spherically symmetric horizon in 
a spherically symmetric spacetime. This would require 
some additional condition to be imposed that selects out 
a unique trapping horizon from the many that occur in 
dynamical black hole spacetimes. 

Despite their quasi-local nature, closed marginally 
trapped surfaces do have some non-local behaviour. 
Their dependence on the slicing is but one manifesta- 
tion of this behaviour. The fact that marginally outer 
trapped surfaces can be found all the way to the event 
horizon [21] is another manifestation. Although the event 
horizon is a fully non-local teleologically defined surface, 
it still acts as the boundary of outer trapped surfaces in a 
Vaidya spacetime that satisfies the null energy condition 
and asymptotes to a static Schwarzschild solution in the 
far future. 

This behaviour is related to both the choice of the 
surface null normals and the requirement that the 
marginally trapped surface be closed. Closed trapped 
surface cannot be found entirely in a fiat spacetime. How- 
ever, parts of a closed marginally trapped surface can be 
found passing through a region of fiat space where both 
null expansions are negative [22]. 

It may be that in the quantum context only an effec- 
tive horizon is meaningful or that the existence of many 
intersecting horizons contribute to the full path integral. 
For purely classical situations one may be able to refine 
the definition of the surface of a black hole. 



V. APPENDIX - SOME USEFUL RELATIONS 
IN THE VAIDYA SPACETIME 

These relations are provided here as a repository 
for reference in the main text. In the advanced null 
Eddington-Finkelstein coordinates (i',r, 0, 0) the metric 



takes the form 

ds^ = -Adi;2 + 2drdv + r^dO^, 



(36) 



where A = 1 - HlIiM A' = ^^, A = -^ and m = 

Radial null vectors, with canonical normalization fn^ — 
— 1 and Ua = —dyr: 



t = (l,y,0,0 



A 



1,0,0 



n" = (0,-1,0,0), 
Ua = (-1,0,0,0). 

Expansions of the null normals: 



A 



Hi) 



Hn) 



(37) 

(38) 

(39) 
(40) 



(41) 



The location of the spherically symmetric MOTS is ob- 
tained by setting 9(^i) = leading to r = 2m{v). The 
variations of the expansions along n" : 



n"V„0(,) . 


A 
r2 " 


A' 
r 


(42) 


^"V„0(,) = 


A 
r 


A2 AA' 
2r2 ' 2r ' 


(43) 


n"V,0(„) = 


2 
= -^ 


1 


(44) 


rV„0(„) = 


A 
= ^■ 




(45) 



Other optical scalars: 

Surface gravities: 

A' 
2 ' 



(46) 



n"rV64 



(47) 

K(„) = -rn^Vtna-O. (48) 

Components of the energy-momentum tensor: 



T„, 



47rr2 



naUb , 



Tnhi 



anh 



47^7-2 
TabU-'t = Tabu-'n^ = . 



(49) 

(50) 
(51) 



Normal to the trapping horizon A = 0: 

7V" = (l,-2m,0,0) , (52) 



iVa = (-2m, 1,0,0) 



(53) 
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Radial tangent to the trapping horizon 

y" = (1,2m, 0,0), (54) 

V; = (2m,l,0,0) . (55) 

Variation of the area along the trapping horizon 

F«V.A = 1^ . (56) 



For this null vector to be the generator of an isolated 
horizon it is required to satisfy TabV°'v'' = giving A = Q 
and therefore from above C = 0. Thus any null vector 
generating an isolated horizon must be proportional to 
Oin) but from the appendix we see 0<n) = — 2/r. This 
only vanishes at infinity. There may be a degenerate 
planar isolated horizon "at" infinity. 



Connection on the normal cotangent bundle 

'I'a = -g>cV6r=(0,0,0,0) . 

Scalar curvature of the horizon 






(57) 



(58) 



The null energy condition is satisfied if £" > 0. This is 
easy to see in a fiducial spherically symmetric null tetrad 
£°',n°', m"', fh°' since for a general null vector w" with 



v" = At + Brf- + Cm" + Cm'', 



(59) 



then Tabv'^v^ = A^Tad^t = ^A^tnaf. The nuh 
energy condition is closely related to whether there is 
positive or negative energy flowing into the black hole. 

If m 7^ there are no isolated horizons in the Vaidya 
spacetime since the norm of the generalized null vector 
(59) is 



V Va 



2AB + 2CC = 0. 



(60) 
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